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Abstract 

The following is a thermodynamic analysis of a III order (and some aspects of 
a IV order) phase transition. Such a transition can occur in a superconductor 
if the normal state is a diamagnet. The equation for a phase boundary in an 
H—T (H is the magnetic field, T, the temperature) plane is derived, by con- 
sidering two possible forms of the gradient energy, it is possible to construct 
a field theory which describes a III or a IV order transition and permits a 
study of thermal fluctuations and inhomogeneous order parameters. 

I. INTRODUCTION 

If the Ehrenfest classification were used to describe a third or a fourth order phase 
transition, the free energy and its low order derivatives, e.g., entropy and the specific heat, 
will be continuous, but, say for a III order transition M, the specific heat will have a 
discontinuous temperature derivative. For a IV order transition, the temperature derivative 
of the specific heat will be continuous, rather, the second derivative of the specific heat will 
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be discontinuous. 



In fact, this is not enough. Consider a II order phase transition |2[] boundary. The phase 



boundary (in the H—T plane, where H is the mechanical variable and T, the temperature) 
is given by 



Here, AC = C Q — Cd where C Q and d are respectively the ordered and disordered phase 
specific heat (C = — TC§) and Ax = Xo — Xd is the similar discontinuity in susceptibility 



transition. A horizontal second order phase boundary may have AC = 0. A higher order 
transition thus necessarily requires both AC = Ax = 0. Stated differently, near the 
transition temperature, F(T) = —{T c — T) n and F(H) = —(H c — H) n . The Ehrenfest 
classification assumes that n = n'. In a scaling description of a II order phase transition, 
n = 2 — a, where a, the specific heat exponent, is small. It is, however, clear that there is 
no limit on what a can be, whether positive or negative. Thus when a is negative and large 
in magnitude, the transition is probably better described as with an order corresponding to 
the nearest integer to n. Any residue in n then can be viewed as a fluctuation contribution. 

The well known liquid-vapor phase boundary for water is mostly I order, except when 
this phase line terminates at the critical point. At the critical point, the transition is second 
order while past the critical point, the transition is continuous without any free energy 
singularities. A possible scenario is that a first order line may terminate at a higher order 
phase transition. We encounter such a possibility in the discussion of magnetic field effects. 

Recent observations in several high T c superconductors have been cited as questioning 
whether the order of the transition is 2. Probably the most extensive study || comes for 
Bao.6K .4Bi03 (BKBO), a cubic superconductor with a T c = 30 K. Here both specific heat Q 
and susceptibility || have been measured and found to be continuous at the transition. There 
are indications HQ that specific heat in Bi and Tl based cuprates has features characteristic 
of a higher order phase transition. The issue still needs to be resolved by careful analysis 
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(X 



fj^r). Having AC = is, naturally, not enough reason to argue for a higher order 



of both the specific heat and the susceptibility. In particular, we need to know whether the 
transition is of order III or IV. 

In the following, we consider both a III order as well as a IV order transition. Section II 
contains a thermodynamic analysis of a higher order phase transition. Section III contains 
a field theory, a Ginzburg Landau type free energy which describes such a transition. This 
section contains only properties associated with a uniform order parameter. A non-uniform 
order parameter is the subject in Section IV, whether caused by an external magnetic field 
or thermal fluctuations. Both are analyzed. Finally, Section V contains a summary and 
some outlook of future. 



II. THERMODYNAMICS 



The objective in this section is to derive an equation equivalent to the Clausius-Clapeyron 
(C-C) equation for a higher order phase boundary and discuss its consequences. There are 
three quantities, ||S , |j| = Q£f and each continuous across the transition. The usual 
derivation of the phase boundary equation proceeds by equating the changes in each of the 
above quantities as one moves along the phase boundary. The resulting equalities are then 
solved for the slope of the phase boundary. Thus, the equations are (making frequent use of 
the Maxwell relation f§ = f^): 

8H Ad 2 S/dT 2 Ad 2 M/dT 2 Ad 2 S/dH 2 Ad 2 M/dT 2 

df ~ Ad 2 S/dHdT ~ ~ Ad 2 M/dTdH ~ Ad 2 M/dH 2 Ad 2 S/dH 2 ^ 

Of all of the possible expressions (some of which can be found in Pippard [[J the one 
below is special for two reasons: Firstly, the numerator contains thermal quantities and 
their derivatives, while the denominator contains mechanical expressions, similar to the C-C 
equation (Eq. (1)) for a second order transition. Secondly, the equation has a symmetric 
form, suggestive of a simple extension to higher order transitions. 

~dH 



Ad 2 S/dT 2 AdC/dT 

Ad 2 M/8H 2 ~ T c Ad X /dH ^ 
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Eq. (3) involves the first derivative of x with respect to H . Thus the phase boundary 
here is determined by the nonlinear susceptibility. If the correction to the susceptibility is 
quadratic in H(x(H) = xo + X2H 2 ) then assuming that specific heat field dependence is 

4 

weak, the phase boundary can be integrated to show that at low fields, AT C oc Ha, a result 
which is different from a conventional superconductor where the Abrikosov result for AT C is 
linear in magnetic field. 

Eq. (3) suggests the form of the phase boundary for a IV order phase transition. One 
might guess the (sometimes called "Ehrenfest equation") phase boundary in a IV order 
phase transition to be 

(dH\ A _ Ad 2 C/dT 2 

[dTj ~ T c Ad 2 x/dH 2 ^ ' 

Going back to III order, given that we expect F = — (T c — T) 3 , the specific heat C(T) is 
given by 

C m (T) = —Td 2 F/dT 2 ~ T(T C - T) (5) 

The specific heat in fact has a broad "peak" at T c /2. A transport measurement might 
report a true T c while a thermodynamic measurement will report a smaller T c . Moreover, 
if, as is customary, the calorimetric measurement of T c is made by deriving T c from a point 
where the entropies (of the order and disordered state) match, one will have yet another 
value for T c . The essential fact is that the transition is not broad, it is described by a 
discontinuity in higher order thermodynamic derivatives. 

Looking back at Eq. (1), we note another curious thermodynamic fact. Since a typical 
metal is paramagnetic and a superconductor diamagnetic, it is not possible to have Ax = 0. 
In fact this condition is satisfied only when the normal state is a diamagnet, as is the case 
with all of the material examples mentioned above. A higher order transitions is possible 
only if the normal state is a diamagnet. 

Finally, we can estimate the thermodynamic critical fields. Thus, in the standard way 

B 2 



2/i 



= (T c - T) n (6) 



i.e., B c ~ (T c — T) n / 2 , a result that is reasonable in the well known II order case but entirely 
unexpected for the higher order phase transitions. In the following we construct specific 
models to better understand the possible microscopic origins of these results. 



III. GINZBURG-LANDAU THEORY 

For a III order phase transition, a Ginzburg-Landau theory contains some surprises even 
though it is relatively straightforward to obtain the characteristic form of the free energy 
Fh(A,T), here A is the order parameter. refers to the free energy with a homogeneous 
order parameter and Fq is the normal state free energy. 

F n = F + aA 4 + 6A 6 (7) 

Here a = ao(|r — 1) and b > a constant. Searching for the minimum of the free energy as 
a function of A, we find 

a9 2\a\ 2an , T. m m 



= T>T C (8) 

The various thermodynamic quantities (all of the quantities refer only to the condensing 
degree of freedom) are given by 

{ > 27 6 2 27 6 2 v Tj [ ' 

° " "9 6 2 T C 2 v 1 T c ) [LL) 

We see that the specific heat has a characteristic temperature dependence. It shouldn't be 
taken for a broad discontinuity of a II order transition since ~ |; where AT is some 
measure of the transition width. 

We should also recall now that there is no quadratic term in the order parameter in 
Eq. (7). In the presence of a quadratic term, the transition becomes a I or a II order ||. 



This is a rather important and subtle point. For example, the magnetic transition in solid 
3 He || involves interaction energies that do not have a quadratic term, they have only a 
quartic term. However, there is always a quadratic contribution from the entropy. Thus 
the transition is usually a I order one. A third order transition requires vanishing quadratic 
term in the free energy. One possible interpretation of Eq. (7) can be that the density of 
states at the Fermi surface (if Eq. (7) is the free energy for a superconductor) which appears 
as an energy scale determining factor, is in fact proportional to A 2 . Thus, we could imagine 
a transition between an insulator and a superconductor. The Fermi surface density of states 
is zero in the normal state indicating an insulator. It is finite in a metal/superconductor 
and the transition would be a curious but a profound feedback phenomena. As far as I am 
aware, there are several quadratic terms in a real insulator-superconductor transition and 
the real transition is a robust II order one. The above analogy is only an illustration. 



What happens if the transition is a IV order? A possible free energy can be written flO 

as 

F h = F + a\A\ 6 + b\A\ 8 (12) 
a = a (T/T c - 1) 

Here, one might imagine that the prefactor density of states at the Fermi surface is driven 
by superconductivity to be quartic in the order parameter. The results corresponding to 
Eqs (8-11) are 

A « = + ^r (1 - T ^ a < < 13) 

W = (14) 



256 b 3 
97 T 

Cw(T) = - ¥f - 2 (l-T/T c r (15) 



IV. SPATIALLY INHOMOGENEOUS ORDER PARAMETER 



Here we consider two effects, first the effect of magnetic field and then the effect of 
thermal fluctuations. To do either, we need to supplement the free energy Eq. (7) with 
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terms that involve a spatially varying order parameter. Again there are many choices and 
it is less obvious how they could be eliminated. The following contains consequences of one 
selection, we will comment on another selection, briefly, at the end of this section. 



A. Model I 



Thus, we supplement the free energy (Eq. (7)) by 

Fq = C|VA| 2 



(16) 



To be specific, we will assume that A(r) is a complex function of position in 3d. Let us 
first look at the temperature dependent correlation length. The Euler-Lagrange equation 
for the order parameter is given by 
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5A 



(F h + Fq) = + 3bA b - CV 2 A = 



(17) 



We look for small amplitude variations of A(r) by writing 

A(r) = A + %) 



(18) 



The linear equation determining <5(r), the final solution being dependent on boundary con- 
ditions, is 



- CV 2 5(r) + { - 



1 5 2 F h 



2 5A 2 



5(r) = 



(19) 



A . 



leading to a length scale, the temperature dependent correlation length £ (T) 



em = c 



[1 S 2 F h 




2 5A 2 


A - 



3bC/S\a\ 



(20) 



Thus, in contrast to the familiar notion of length scale, here there is an asymmetry about 
T c . For T > T c , the correlations are non-exponential. Below T c , the correlations length 
depends linearly on l/|a|, as opposed to the square root dependence in a second order phase 
transition. 
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We can also identify the superfluid density. A current can be viewed as that arising from 
the space gradients. In that case, the density is given by expressing the free energy due to 
current as VzPsVg where v s is the supercurrent. With Eq. (16), we have p s ~ A 2 and the 
consequent temperature dependence. 

A.l. Magnetic Field Dependence 

The gauge invariant free energy requires that the charge interaction with magnetic field 
be described by a transformation of the gradient term V — > (V + ^^4) here 0o is the flux 
quantum h/2e and A is the vector potential (B = V x A refers to the local field, which 
approaches the external field H outside). 

It is straightforward to show that the equation for the vector potential represents flux 
expulsion whenever A 7^ 0. It is also true that in analogy with a conventional supercon- 
ductor, \~ 2 oc A 2 . All these features are standard and are dependent on n s = A 2 . The 
difference appears, not entirely unexpectedly, at the calculation of H C 2(T). The curious 
result is that the transition becomes I order in the presence of a magnetic field. 

To see that, we derive the criterion for the instability of the normal state. The lowest 
order Euler-Lagrange equation is given simply by 



The solutions of this equation are dependent on the choice of gauge, but the eigenvalues 
A = (2n + l)^-H is independent of gauge choice, with the lowest energy level corresponding 
to n — 0. If the corresponding wave function is 0„(r), we can write 



where n stands for the set of quantum numbers necessary to characterize the state. In the 
mean field sense, the free energy can be written as 




(21) 



Afe) = ECn0n(r) 



(22) 



(n) 



2nC 



HQ + aCo 4 + &Co 



(23) 



F = 



0o 



The transition occurs at H c2 (T) where 
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HAT) = ^ (24) 

However, the order parameter develops discontinuously at T c . Its value at T c , A( is 
given by 

M^H* < 25 > 

The latent heat (L = T c As) can also be obtained, it is given by 

L = T c As = 2 -^H/T c (26) 

All of these quantities are measurable and provide an unequivocal test of the model for the 
gradient energy. They represent a first order phase transition for H ^ while the H = 
transition is a III order one. We have a first order line ending in a critical point which is of 
higher order than the conventional second order. 

At this stage, the applicability of ideas of this subsection to BKBO breaks down. As far 
as we can tell, the transition is not of I order in any finite field, there is no latent heat, no 
discontinuity — nor any other trace. In fact, all evidence in finite field, points to a transition 
higher in order than II. Before proceeding with another ansatz for free energy in sec. 4.2, 
let us consider thermal fluctuations following Eq. (16). 

A. 2. Thermal Fluctuations 

Let us consider a scalar order parameter A(r). We consider a partition function 

z = J D[A(r)]exp[—P(F h + F G )] (27) 

The partition function provides the thermodynamic free energy which then leads to the 
other thermodynamic properties. As is customary, the integral is evaluated in a saddle 
point approximation. Consider A(r) = A + 5(f) where A is the minimum of F h , as in 
Eq. (8): 



e-P f ' = z = e-W>l[Jd6 k exp 



(3 y, f S 2 F h 



2 V U A2 



+ ^ 2 61 

A > 



(28) 



where 5 q (or 5 k) is the Fourier transform or 5{r), and the first term is given by 8|a| 2 /36. 
Thus, the free energy becomes 



F = (Fh) + F f 
Ff = - k ^n 
5 2 F 



(29) 



C, = -T^ - \a\ d ' 2 (30) 

In this model, the transition at H = is quite robust. In Eq. (30), we have the temperature 
dependence of the fluctuation specific heat, which is well behaved for d > 2. The upper 
critical dimension for this model is d = 2. (Recall that the upper critical dimension in a II 
order case is d = 4). The fluctuations in this model, described by Eq. (16) are divergent 
only for d — 1. Thus there are no critical fluctuations and no need to go beyond mean field 
theory, at least not for d > 2. This also means that if the experimental exponents are not 
meanfield like, model I is probably not applicable. 

If A(r) were a 2-d vector, then the phase fluctuations behave similar to their behavior 
in case of a II order phase transition. There are the usual infrared divergences in 1 and 2-d 



indicating the importance of topological defects [11] in the phase transitions, etc 



B. Model II 

Let us consider a different gradient energy term. If we follow the physical picture where 
the overall energy scale in the free energy depends on the order parameter, the gradient 
term perhaps looks like 

F G = 5A 2 |VA| 2 (31) 

which can also be viewed as the square of the gradient of A 2 . We see that for a scalar order 
parameter, the Eyler-Lagrange equation is given by 

- 2|a|A 2 + 36A 4 - g(VA) 2 - £AV 2 A = (32) 
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A linearized version of this equation, about the homogeneous solution, becomes 

A\a\8(r) - 5V 2 5(r) = (33) 

In contrast to Eq. (19), the correlations length exponent is back to being \. But much else 
has changed. The effect of magnetic field can be incorporated by replacing the derivative 
in Eq. (31) by a covariant term. Like sec. IV. A. 1, the order parameter can be expanded 
in the Landau orbital wave functions which are eigen functions of the covariant Laplacian. 
The nonquadratic nature of the free energy introduces interactions and transitions between 
the Landau orbital states. However overlooking these interactions (for small fields as in 
sec. IV.A.l), we see that the transition in the presence of a magnetic field is III order with 
a phase boundary given by 

lire 

-^H = a (l-T/T c ) (34) 

If A(r) is a complex number, the number density of excitations is proportional to Aq, a 
feature also shared by the penetration depth (A -2 oc A 4 ). 

The proverbial fly in this ointment comes from a study of thermal fluctuations. The 
thermal fluctuations contain a divergent (logarithmic) contribution to the free energy, yield- 
ing a specific heat that diverges as |a|~ 2 , independent of dimension. There is additional 
dimension dependent divergence, which for specific heat takes the form C oc | « | C^/ 2 2 ) _ To 
see the details, consider Eq. (28) as worked out for model II, we have 



e -f>fi = z = e -P(FG)Y[ [ d5 k exp 



-/3£A 2 [4|a| + cg 2 K 



(35) 



kT 

i.e. F f = — 5>[A 2 (4|a|+cg 2 )] ( 36 ) 
<? 

It is the Aq term that is singular in the eventual derivation of the entropy and specific heat. 

At this moment, it is hard to imagine that the entire mean field analysis of model II, as 
described above, is meaningless. And yet, a mean field analysis depends on the validity of 
the saddle point evaluation of the partition function. An estimate of corrections to mean 
field, seen above for model II, is divergent and it is unclear if some sort of renormalization 
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will restore these results to a finite value. Barring that, model II remains very dubious as a 
valid model. These problems are equally exacerbated for a IV order phase transition. 

V. SUMMARY 

It appears therefore that a thermodynamic description of a III or IV order phase transi- 
tion is relatively straightforward, although the consequences are almost unavoidably subtle. 
Once we recognize the possibility of these phase transitions, it is clear that we cannot lump 
them together with all of the other continuous transitions. There are subtle and interesting 
differences which call for a separate identification and analysis. 

It is possible that a higher order transition has been observed before. Then it was 
attributed to an extreme case of sample inhomogeneity since the transition temperature 
measured from different techniques were all different. The above is an attempt to outline 
a systematic formalism for the analysis of higher order transitions. We have an equation 
for the phase boundary Eqs. (3) and (4) and furthermore, we also have a field theory which 
leads to a III (or IV) order phase transition. In searching for the effects of a magnetic field or 
thermal fluctuations, we find two models. In one case, the finite field transition is a I order 
transition which ends into a III order critical point. All thermodynamic quantities vanish 
at this critical point (H = 0) appropriately. In the other model, we find strong thermal 
fluctuations which, on the one hand appear to question the validity of a mean field theory. 
It is also possible that the correct model still remains elusive. 
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